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ABSTRACT 

We report on the detection of reddening toward z ^ 0.2 galaxy clusters. This is measured by correlating 
the Sloan Digital Sky Survey cluster and quasar catalogs and by comparing the photometric and spectroscopic 
properties of quasars behind the clusters to those in the field. We find mean E{B — V) values of a few x 1 0^^ mag 
for sight lines passing ^Mpc from the clusters' center The reddening curve is typical of dust but cannot be used 
to distinguish between different dust types. The radial dependence of the extinction is shallow near the cluster 
center suggesting that most of the detected dust lies at the outskirts of the clusters. Gravitational magnification 
of background z ^ L7 sources seen on Mpc (projected) scales around the clusters is found to be of order a few 
per cent, in qualitative agreement with theoretical predictions. Contamination by different spectral properties 
of the lensed quasar population is unlikely but cannot be excluded. 

Subject headings: dust, extinction — galaxies; clusters: general — intergalactic medium — quasars: general 



1. INTRODUCTION 

Zwicky (1957) suggested the presence of dust in the intr- 
acluster medium (ICM) to explain galaxy number counts be- 
hind the Coma cluster and a search has begun to quantify its 
properties. Similar analyses were carried out by Karachentsev 
& Lipovetskii (1969), Bogart & Wagoner (1973), and Boyle 
et al. (1988 using quasars) and yielded B-band extinctions 
of order 0.2 mag. However, a more recent study by Maoz 
(1993) using radio-selected quasars behind rich Abell clus- 
ters yielded only upper-limits on the reddening [E{B — V) < 
0.05 mag] with similar limits [E{B — V) < 0.02 mag] obtained 
from the analysis of galaxy colors behind APM clusters (e.g., 
Nollenberg et al. 2003 see also Ferguson 1993). In contrast, 
Hu (1992) found E{B — V) ~ 0.2mag toward cooling flow 
clusters (but see Annis & Jewitt 1993). Dust emission from 
clusters has been marginally detected in the infrared resulting 
in a dust-to-gas ratio as low as ^ 10^^ in a few Abell clusters 
(Stickel et al. 2002; but see Bai et al. 2007). 

Cluster members are known to contain dust (e.g.. Edge et 
al. 1999, Popescu et al. 2002, Boselli et al. 2004) and galax- 
ies are seen to shed gas and dust into their environments (e.g., 
Heckman et al. 2000; see also Jonsson et al. 2006). However, 
the survival time for small dust grains mixed with the ICM 
is much shorter than Hubble time (Draine & Salpeter 1979), 
and would be absent from the cluster core (but could survive 
on large scales; see Dwek et al. 1990) unless constantly re- 
plenished (e.g., Montier & Giard 2004) or shielded from the 
hot plasma (e.g., Voit & Donahue 1995). Presently, there are 
few meaningful constraints on the composition and quantity 
of intracluster dust, an understanding of which would have 
important implications for cluster evolution and metal enrich- 
ment processes, as well as for dust physics. In this Letter 
we provide new answers to 50 years old questions concerning 
dust in galaxy clusters, by harnessing the power of the Sloan 
Digital Sky Survey (SDSS). 
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2. METHOD 

We use a 0.1 < z < 0.3 sample of ^ 10^ optically-selected 
clusters from Koester et al. (2007) with velocity dispersions 
>400kms-i. The mean /?200 (the radius within which the 
galaxy density is 200 x the field density) is ~1 Mpc^. The 
cluster center is defined by its brightest member whose posi- 
tion may deviate from the center of mass by a fewx 100 kpc. 
This sample was correlated with the SDSS/DR5 spectroscopic 
quasar catalog (Schneider et al. 2007), identifying all quasars 
within 10 Mpc of the clusters' centers. All data were cor- 
rected for galactic extinction using the values of Schneider 
et al. (2007). In particular, the spectra were corrected using 
the reddening curve of Cardelli et al. (1989). 

We employ a statistical approach to search for dust by 
comparing the photometric and spectroscopic properties of 
quasars at different impact parameters, b, from the cluster 
center with those of quasars in the field. The SDSS quasar 
selection criteria have no effect on the results and do not bias 
the observed colors for E{B — V) ^ Q.lE{g — i) < 0.1 mag 
(see Menard et al. 2007). As clusters vary in size, we define a 
normalized impact parameter: 



b = b/R 



200- 



(1) 



A cluster would induce a mean color excess in background 
objects whose sight lines pass a distance b from its center (see 
appendix A): 

-'^^'^ = ^^''- + ^^^i + (2) 

where i,j denote SDSS photometric bands, T the ef- 
fective dust optical depth, and 5^ the foreground 
emission by the cluster These terms result in effec- 
tive reddening each having very different signatures 
(see Fig. 1 where median quasar/early-type cluster 
galaxy colors were used so that [M,^,r,/,z]quasaiy galaxy = 
[19.4, 19.2, 19.1, 19.0, 18.9]/[21.2, 19.9, 18.8, 18.3, 18.0]). 
The ^(/x) term arises when quasar colors are luminosity 
dependent and gravitational magnification, pL is considerable 
(see §3 and the appendix). 



We use the standard cosmology with Hq ■ 
and = 0.7 



^70kms-'Mpc-', il,,, =0.3, 
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Fig. 1 . — (^(g — i)) as a function of b is consistent with a h ' decline at 
large radii and a flat core at smaller radii (see overlaid dashed line curves). 
The inset shows the results for a particular annulus (b ~ 5) for all SDSS 
colors, u — g is considerably noisier than the other colors (see text). The 
color excess for g — r, r — i, / — z is consistent with known reddening curves 
(solid line) indicating that foreground emission (dashed line) is negligible 
(see also §4 and appendix A). 

The analysis proceeds as follows: we define an annulus 
about the cluster center in the range [b,b + 5b]. All quasars 
within this annulus serve as a primary sample (consisting of 
^ 3000 objects for b ^ [0, 1] with the number oc b5b). Many 
primary samples span the range b G [0,7] with 5b being set 
by the requirement for adequate S/N and spatial resolution. 
We define our control sample to be all quasars with b>7. As 
the control sample is large, we choose many random control 
sub-samples with the same number of objects as the primary 
sample and compare the geometric mean of their spectral and 
photometric properties. This allows us to quantify the mean 
color excess and its variance via Monte Carlo simulations. 

For the photometric analysis we choose individual control 
sample quasars to have the same redshift (up to ±0.1) as that 
of individual primary sample quasars. This ensures that we 
are not sampling slightly different observed redshift distribu- 
tions caused by gravitational lensing. Similarly, for the spec- 
troscopic sample we work in narrow redshift bins {5z — 0.2) 
so that our results are not masked by the smearing of quasar 
emission lines to continuum-like features (e.g., York et al. 
2006). To test for potential biases we verified that both anal- 
yses produce a null mean color excess and a mean spectral 
ratio of unity, respectively, when comparing different control 
sub-samples. 

3. RESULTS 

The color excess is shown in Figure [T| as a function of b. 
Clearly, reddening in g — i is detected with high significance. 
The color shifts are consistent with known reddening curves 
so that {E{g-r)) > {E{r-z)) > {E{i-z)), indicating that 
foreground emission by the cluster is negligible (see Fig. [T] 
and §4). Our {E{u~g)) estimate is less secure given the 
known calibration issues associated with the M-band filter^. 
We find (E{g — /)) values of a few x 10^^ mag. Inspection of 
the color-excess distributions (not shown) reveals that the re- 
sults are not dominated by a few outliers and that most quasars 
are redder behind clusters. The qualitative agreement with 
known reddening curves suggests that £, is negligible. 

Motivated by our photometric results we assume that the 
flux of background quasars is subject only to extinction and 
magnification so that the ratio of the observed (flux) spectrum, 
Fobs- to the intrinsic one, Fi„t, is 

FoU^)/FUX)=^ie-''-^^ (3) 
* see e.g., |http://www.sdss.org/dr5/algorithms/fluxcal.html| 



where t(A) is given by the extinction curve. In our spectro- 
scopic analysis we consider 1 .5 < z < 2 quasars which ensures 
i) that lensing effects, hence ^ , are smaller (see appendix), ii) 
adequate S/N at each redshift bin (see Fig. 3 in Schneider et 
al. 2007 showing that ~ 27% of DR5 quasars are in this red- 
shift range), and iii) adequate wavelength coverage excluding 
the Balmer continuum and emission lines band (which we ver- 
ified contaminate the signal, especially in lower-z objects; Yip 
et al. 2004). A representative case is shown in Fig. |2] which 
is nicely fit by a Cardelli extinction curve (other curves give 
similar results). Deviations from a smooth reddening curve at 
the location of the C IVA1549 line are probably a manifesta- 
tion of the Baldwin (1977) effect (see appendix B). Using 
statistics while avoiding emission line regions and the Balmer 
and Fe II band, gives us the reddening and the extinction, al- 
lowing us to deduce the extinction-corrected flux ratio i.e., the 
mean observed gravitational magnification, /Xobs (see Fig. |2]). 
We carried out the analysis for each redshift bin and averaged 
our results over all bins to produce Fig. [3] Typical {E{g — /)) 
values are larger than the photometrically deduced values on 
account of the latter being contaminated by the Balmer/ Fe II 
band. The observed magnification is of order a few per cent. 

4. CONCLUSIONS 

We deduce finite extinction for quasars behind clusters and 
our results are consistent with a large covering factor of dust 
toward these objects. The dust mass enclosed within a sphere 
of radius b is 

M,..^10^m'f|^M,. (4) 
\Mpc / 10 ^mag 

(Here we used the ICramers-Kronig relations and grain density 
of 2.5g cm^^; Kxiigel et al. 2003). This implies a dust-to-gas 
ratio for the ICM which is < 5% of local inter-stellar medium 
(ISM; here we used ICM mass of 4 x lO'^M© within the cen- 
tral Mpc; Morandi et al. 2007). If, as suggested by some 
theoretical models, galaxies eject their ISM in amounts com- 
parable to their own mass during their formation, then our re- 
sults imply that either the associated dust is in a clumpy form 
with a small covering factor or else it is rapidly destroyed and 
the metals are deposited in the ICM. Large scale distributions 
of dust could arise from a non-virialized gas component at the 
outskirts of clusters left from past merger events or related to 
halos of neighboring galaxies and small groups provided their 
size is of Mpc scales. Singular isothermal sphere (SIS) mod- 
els, where dust follows dark matter, seem to over-predict the 
extinction on small scales (see Figs. [U&O indicating that the 
centers of clusters may be relatively devoid of ISM-like dust. 

Our results indicate that foreground emission is not impor- 
tant. This means that the contribution of foreground light that 
enters a (median) quasar's point-spread function is > 7 mags 
fainter in the z-band (Fig.l and appendix A) and requires that 
most sight lines pass > 10 (de-Vaucouleur) half-light radii 
(~ 15kpc; see Fig. 4 in Bernardi et al. 2003) from a me- 
dian cluster member. Repeating the entire analysis for quasars 
whose sight lines pass > 30kpc from confirmed members 
yields similar results. We estimate foreground light contam- 
ination by resolved and non-resolved cluster members (using 
a luminosity function with a slope of —0.8 extrapolated down 
to ^ 0.04L*, and including 25% of blue galaxies; Hansen et 
al. 2007) to be at the per cent level, and hence negligible. 

The observed magnification was found to be few per cent 
on Mpc scales. Correcting for luminosity function effects (see 
appendix B) the intrinsic magnification is of order 10%, which 
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Fig. 2. — The mean spectral ratio for quasars with 1.7 < z < 1.9 and at various annuli around the cluster center to field quasars (h > 7) at the same redshift 
bin (the primary samples with b G [1.5,3], [3,6] consist of 2 x 10^, 8 x lO' objects, respectively). In all cases the ratio is well fitted by a reddening curve (green 
line). The dashed red line shows the flux ratio after correcting for the effects of extinction and approaches unity for large h, as expected (see text). Our algorithm 
does not introduce an artificial signal as seen in the right panel (a low-level signal is observed since the control sample is composed mainly of > 8 quasars). 
Prominent quasar emission lines are marked as is the wavelength beyond which the Balmer continuum and Fe II line bands are present (see text). 



Fig. 3. — Reddening and magnification from the spectroscopic analysis as 
a function of scale b. Left panel: Positive E(g — i) (i.e., reddening) is detected 
with high significance. Here too, reddening rises less sharply toward the cen- 
ter compared to a 5^ ' law and is consistent with a flat core profile (overlaid 
dashed line curves). Right panel: The observed magnification, /Id,, — 1 as a 
function of b is of order a several per-cent. No significant deviations can be 
seen from a magnification profile (dashed line). 

is in qualitative agreement with SIS model predictions for the 
magnification of a z = 1.7 quasar behind a rich z — 0.2 clus- 
ter being ^ - 1 0.2(/7/Mpc)"' (e.g., Menard et al. 2007). 
Given the uncertainties on the dark-matter profile, our approx- 
imation for the SDSS selection function and the luminosity 



function of quasars, we find this agreement to be surprisingly 
good. It also implies that grey opacity (by large grains that 
survived for a Hubble time at the cluster center) is not re- 
quired. The above agreement, the fact that the deduced red- 
dening and magnification follow different trends at small b 
(Figs. [T|&[3]l, and the good fit to known reddening curves, all 
suggest that our results are not significantly contaminated by 
luminosity-dependent intrinsic color differences in the lensed 
quasar population. The agreement between the photometric 
and spectroscopic results, which are prone to different sys- 
tematics, is also reassuring. 
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APPENDIX 

A. COLOR VARIATIONS IN OBJECTS BEHIND GALAXY CLUSTERS 

Consider how an intervening cluster affects a background quasar's flux, F so that the measured flux 

F' =F^e-^ +Ffl« m,=m, + j^T,-2.51og(^) + 5f (Al) 

where m (m) is the (measured) observed magnitude, jj, is the gravitational magnification, T the (dust) optical depth through the 
cluster, and F-^^^ the contribution from a foreground source associated with the cluster. Here, / denotes the relevant SDSS filter 
band and 5/^ (R) the contribution from foreground cluster emission where 

ff/g , ff/g 
5f^^M-V^^. (A2) 

Defining the mean color of quasars at a given impact parameter b as {i — i)^, = {mi{b) —mj(b)'j, then their color excess with 
respect to field quasars (i.e., atb ^ °°) is 

(£(/ - ^ (/ - j), {i ./■)„ = ^{Ti- T/), + {5t- 5f ), + (A3) 
The color excess by foreground emission may be recast in the form 



(5f -5/) =C(m,-,mj)/"'' where C = 2.5 



and Am; = ( m{^^ -nii) . (A4) 



Here m^J^ is due to foreground emission (see the inset of Fig. [T] where a Am^ — 7 model is plotted in dashed line and a typical 
early-type galaxy spectrum is assumed to characterize the cluster's spectral energy distribution as a whole; see §2). The last term 
in equation IA3 I is: t,y = (m, — Wj)^ — ('w, — Wj)^ which stands for the "intrinsic" color difference between quasars behind the 
cluster and field quasars. This effect can be considerable (/(a) gravitational lensing is important and (b) quasar colors depend on 
their intrinsic luminosity. The origin of this effect is detailed below. 

B. MEAN PROPERTIES IN MAGNITUDE-LIMITED SAMPLES 

Consider a property Q of quasars which depends on their intrinsic luminosity (or magnitude; assume for simplicity a narrow 
redshift range). In a magnitude limited sample with selection function 0(m — m/) (here we approximate to be a step function 
with a limiting magnitude m/ = 19.1; Schneider et al. 2007; see also Menard 2005), the average of Q is 

{Q)=C^^J dm@{m-mi)Q{myV{m) ; C^v = j dm@{m-mifV{m). (Bl) 

Here ^ is the luminosity function. The difference between the measured {Q) when magnification operates (i.e., m m + 
5m; dm < 0) to when it is not, is to first order 

5(e)~C^>5m r dm&{m~ni,)(^^> + Q^] ^ ^ ^ ~ (a + j3)5m (B2) 

J-oo \dm dm J {Q) 

where in the last step we considered a specific example in which Q{m) ^ e""\ and ^'{m) °^ e^"\ Clearly, the effect on mean quan- 
tities can have either sign, depending on the sign of a + j3. Taking, for example, a — —0.8 which characterizes the dependence 
of quasar emission line luminosity, L/, on the quasar luminosity (via the Baldwin relation; e.g., Dietrich et al. 2002) and /3 = 1 
for the quasar luminosity function, we find that 5 (L/) / (L/) ~ — 1% which is in qualitative agreement with our findings (note the 
feature at the C IV A 1549 wavelength; Fig. |2]i. Less elegant expressions for geometric means lead to similar results in this case. 
Given the uncertainties on a, j3, it would not be surprising to find a different trend for other lines at other redshifts and luminosity 
ranges. Similar considerations apply also for the Balmer continuum and Fe II lines band yet the luminosity dependence is not 
well known. As for measuring the magnification itself, this corresponds to a = — 1 which, for j3 = 1, yields 5 (Q) = (i.e., 

jj, = 1). Clearly, this is not the case and a more realistic luminosity function to consider deviates from a powerlaw and is of the 

-1 



form ^'{m) < 



jQ-/3(m-m*) _|_ jQ-y(m-« 



with = 19.1, j3 = 0.98, and 7= 0.15 (e.g., Myers et al. 2003). Using that, we 



obtain that the observed mean magnification is related to the real one by — 1 ~ 3 x (^obs — 1 ) (see also Menard 2005). 



